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We have explored the spatial distribution of an ion cloud trapped in a linear octopole radio- 
frequency (rf) ion trap. The two-dimensional distribution of the column density of stored AgJ was 
measured via photofragment-ion yields as a function of the position of the incident laser beam over 
the transverse cross section of the trap. The profile of the ion distribution was found to be dependent 
on the number of loaded ions. Under high ion-loading conditions with a significant space-charge 
effect, ions form a ring profile with a maximum at the outer region of the trap, whereas they are 
localized near the center axis region at low loading of the ions. These results are explained quan- 
titatively by a model calculation based on equihbrium between the space-charge-induced potential 
and the effective potential of the multipole rf field. The maximum adiabaticity parameter Tjmax 
is estimated to be about 0.13 for the high ion-density condition in the present octopole ion trap, 
which is lower than typical values reported for low ion densities; this is probably due to additional 
instability caused by the space charge. 

PACS numbers: 37.10.Ty, 41.90.-|-e, 36.40.Mr 



I. INTRODUCTION 



Linear multipole radio-frequency (rf) ion traps have 
become increasingly important in atomic, molecular, and 
cluster physics [ll, l2i- In spectroscopic studies, for ex- 
ample, these traps improve the resolution by enabling 
cooling with a neutral buffer gas and improve the sensi- 
tivity by enabling accumulation and thereby increasing 
the number density |3l-ll4| . In gas-phase reactivity stud- 
ies, multipole rf traps have proved ideal in thermalizing 
all degrees of freedom of reagents and in providing ac- 
curate knowledge of their concentrations [15l - (l7| . The 
success of these traps is due to the characteristics of the 
time-averaged potential experienced by the trapped ions, 
the effective potential 'J'cff- The strength of the effective 
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potential is proportional to p^ for a 2p-pole trap, and 
its dependence on the distance r from the central axis 
is proportional to r^P~^ [18[. For high values of p, traps 
are particularly deep and have a large nearly field-free 
region around the axis, which guarantees little rf heat- 
ing [l8l - l2ll |. The deeper a trap, i.e. the larger the phase 
space acceptance, the larger the fraction of ions that can 
be captured and accumulated. Additionally, buffer gas 
cooling allows phase space compression, which, in turn, 
can further increase the number of stored ions. Recently, 
the depth of a 22-pole trap was measured by analyzing 
the evaporation rates of trapped ions from a thermody- 
namical point of view; this result was employed to calcu- 
late effective trap depths for several 2p-pole ion traps as 
a function of the rf amplitude [22., 23] . 

Taking advantage of the high ion density attainable in 
a multipole ion trap, we have recently been able to ap- 
ply cavity ring-down spectroscopy directly to the mass- 
selected ions stored in an octopole trap [j-lal- In the 
course of these studies, the ions were found to be dis- 
tributed not uniformly inside the trap, as it was first 
pointed out in Rcf. [3]. In fact, the ion density distribu- 
tion is not defined by the effective potential alone in a 
regime where space charge effects become relevant. Un- 
derstanding these distributions in linear rf traps is ex- 
tremely important (i) for the optimization of spectro- 
scopic methods that profit from the maximal overlap of 



laser light with the ion clouds inside the trap, (ii) for 
the measurement of absolute absorption cross sections, 
(iii) to accurately estimate the thermalization tempera- 
ture in reactivity studies, and (iv) for optimization and 
interpretation of experiments where the ion trap is used 
as a pick-up cell [24j . 

The density distribution of a non-neutral plasma in a 
Pennin g tr ap has been extensively studied by Dubin and 
O'Neil [23- For a Paul trap, i.e., a quadrupole rf field, 
profiles of an ion cloud were studied for atomic metal 
ions and organic molecular ions by monitoring fluores- 
cence and photodissociation of stored ions, respectively 
|26l - l28l |. These profiles have a Gaussian shape with the 
maximum ion density at the center of the trap. More 
recently, even single ions have been observed by fluores- 
cence imaging of Coulomb crystals formed in a linear 
Paul trap at temperatures around 10 mK [23]. For mul- 
tipole ion traps, on the other hand, only a few measure- 
ments have been reported. Walz et al. measured a radial 
distribution of ions in a three-dimensional hyperbolic oc- 
topole ion trap by monitoring fluorescence intensity from 
stored Ba"*" ions [30| . They found that Coulomb repulsion 
between stored ions resulted in two separate ion peaks. 
Wester and coworkers reported a radial distribution of 
column densities of 0H~ ions in a 22-pole ion trap from 
photodetachment-rate measurements J3l|. The distribu- 
tion showed a rather uniform profile and was explained 
by a model calculation neglecting Coulomb-repulsion ef- 
fects; this model was applicable to the measurement per- 
formed at a low ion density (less than 10^ stored ions). 
They extended the measurement to two-dimensional to- 
mography of the column density [32| . A Coulomb crystal 
of laser-cooled Ca"*" ions has recently been formed in a 
linear octopole ion trap by Okada et al.; the ions were 
observed by fiuorescence imaging for storage of up to 10* 
ions, which is still in the low-density regime |33l. |34| . 

In this paper, we report on the measured ion-density 
profile in a linear octopole ion trap at a high-density 
regime, where the space charge plays a significant role. 
Silver dimer cations, Ag^, are detected via photofrag- 
mentation and the two-dimensional distributions of col- 
umn densities are measured as a function of the number 
of loaded ions and as a function of the amplitude of the 
rf field. Up to about 10^ ions are loaded into the trap, 
which is the space-charge-limit condition for the present 
ion trap geometry. The distribution profiles are com- 
pared with a model calculation, which takes the balance 
between the trapping force due to the rf field and the 
Coulomb repulsion among the stored ions into account. 



II. EXPERIMENTAL PROCEDURES 



cluster-ion source. A typical current of the Agj ion beam 
after the mass filter was 200 pA, with QMl operated at 
relatively low mass resolution (m/Am ^ 20). The se- 
lected AgJ ions were transported by octopole ion guides 
(IGs) and quadrupole deflectors (QDs) to a linear oc- 
topole ion trap. 

The ion trap has a length of 40 cm and an inner di- 
ameter of 1.1 cm. The trap is fllled with He buffer gas 
to thermalize all degrees of freedom of the ions. All data 
presented here were measured at room temperature. The 
rf potentials at about 3 MHz applied to the octopole's 
rods were provided by a homemade rf generator [35|. 
The rf amplitude Vri was varied in the range 95-300 V, 
while the offset dc potential was kept at —12 V. The ions 
were conflned in the longitudinal direction using poten- 
tials generated by an entrance (Ljn) and an exit (Lout) 
electrode. The potential of Li„ was switched between —8 
and -|-5 V for loading and trapping, respectively. Lout 
was held at -1-5 V during loading and trapping, and was 
switched to —12 V for extraction of trapped ions. The 
number Nq of ions stored in the trap was measured via 
monitoring the ion current at the detector (IDl), taking 
into account the transmission probabilities of IG and QD. 
These were estimated to be unity and 0.9, respectively. 
By controlling the loading time, the ion current, and the 
rf amplitude, Nq was varied from 4.0 x lO'' up to 1.2 x 10^; 
the latter is the maximum ion loading permitted by the 
present ion trap. 

The column density of the ions was measured by 
recording the photofragmentation yield introduced by 
UV laser pulses at 415 nm, near the peak absorption 
of Ag2 |36[. These were generated by using an opti- 
cal parametric oscillator (MOPO-HF, Spectra Physics) 
operated at a repetition rate of 10 Hz, which was attcn- 
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A schematic diagram of the experimental setup is 
shown m Fig. [U a detailed description can be found else- 
where [33. Ag^ ions were selected by a quadrupole mass 
fiher (labeled QMl, model MAX-4000 by Extrel CMS) 
from the distribution generated by a magnetron-sputter 



FIG. 1: A schematic diagram of the experimental setup. QMl 
and QM2, quadrupole mass filters; IGs, octopole ion guides; 
QDs, quadrupole deflectors; FDin and PDout, photodiodes; 
PM: a power meter; Lin and Lout, entrance and exit electrodes 
of the ion trap; IDl and ID2, ion current detectors. 



uated to about 20 ^J/pulse. With a collimated beam of 
2 mm diameter, the intensity was sufficiently low that 
the photofragment yield depended linearly on the laser 
pulse energy. The laser beam was aligned carefully to be 
parallel to the axis of the ion trap. To map the ion den- 
sity distribution, the laser position was displaced both 
vertically and horizontally with an interval of 0.5 mm. 
The laser pulse intensity was monitored during the mea- 
surement by the photodiodes (PDjn and PDout) located 
before and after the ion trap, outside the vacuum cham- 
ber. The signal intensities of the PDs were calibrated by 
a power meter (PM) placed downstream. The number of 
photons of the injected laser pulses, Np, was determined 
from the signal intensities of PDin measured shot by shot. 
The intensity ratio between the two PDs was used to con- 
firm that the laser beam was not clipped while scanning 
its position. The number of photofragment ions, Nf, are 
extracted, mass analyzed using a quadrupole mass filter 
(QM2) and detected at ID2. The transmission probabil- 
ity of QM2 was estimated to be 0.4. 

The measurement was performed following this proce- 
dure: First, AgJ ions were loaded into the ion trap for 
a duration between 0.3 and 2.0 s depending on the num- 
ber of ions to be stored. Second, the stored ions were 
thermalized by collisions with the He buffer gas at room 
temperature for 0.5 s. Third, the ions were irradiated 
with the laser pulse for 1 s, i.e., 10 shots. Finally, the 
ions were extracted from the trap and the yield of Ag+ 
photofragments was recorded. The above measurement 
was repeated five times at each laser position. 

An absolute value of the local number density of ions, 
n{r), was derived from the relationship: 



n{r) 



fir) 
aL ' 



(1) 



where a is the photodissociation cross section, L is the 
length of the ion trap, and f{r) is the number of frag- 
ments ions per photon (Nf/Np). We measured the quan- 
tity / as a function of the laser position r. Because the 
trap is 40 cm long, we assume a uniform distribution 
of ions along the longitudinal direction of the linear ion 
trap. The cross section a was determined from the nor- 
malization condition: 



Nn = 



L / n{r)dr = — f{r)dr. 



(2) 



III. RESULTS 

Figure [2] shows the result of a two-dimensional scan 
of an ion cloud containing 5.5 x 10^ AgJ ions. A major 
part of the trapped ions are found in the outer region 
rather than at the center of the ion trap. The small vari- 
ation of the densities depending on the azimuthal angle 
might be caused by imperfect configuration of the poles 
and/or other neighboring electrodes [32] . In the following 
discussion, we assume cylindrically symmetric distribu- 
tions and will analyze the radial distributions obtained by 



one-dimensional scans along the horizontal axis. These 
data are then analyzed using an adiabatic approxima- 
tion, which provides a rotationally symmetric shape for 
the effective potential p^ . 

The radial distributions of the ion density are shown 
in Fig. |3] for three different amounts of loaded ions, along 
with model calculations discussed in the following sec- 
tion. For these measurements Vrf was held constant at 
200 V, the optimal value. With Nq = 4.0 x 10^, the 
ion distribution is found to be concentrated around the 
center of the trap (Fig. [3^a)). As A^o is increased, the 
ion density in the central region increases only slightly 
while most of the ions are found in the outer region. 
At the maximum loading condition of Nq = 1.2 x 10^, 
the ion density is peaked around r ~ A mm (Fig. EJc)). 
The fact that the density distribution is confined to the 
center of the trap when the ion number is low suggests 
that the ions are well thermalized by collisions with the 
He buffer gas. Therefore, we interpret the ring profile 
of the distribution measured for the highest ion density, 
Nq = 1.2 X 10^, as the effect of Coulomb repulsion be- 
tween the trapped ions. 

The radial distribution profiles are shown in Fig. 2] for 
three values of the rf amplitude K-f. Ions were loaded 
until saturation for each V^f, the number of ions was 
measured to be 1.4 x 10^, 1.2 x 10^, and 1.1 x 10^ for 
Vrt = 95, 200, and 300 V, respectively. The number of 
ions increased by a factor of 8.6 when V^f was changed 
from 95 to 200 V, while it slightly decreased when V^f was 
changed from 200 to 300 V. Although the total number 
of ions varied only slightly between V^f = 200 and 300 V, 
the distributions exhibited clearly different profiles; with 
the higher rf amplitude the density peak became sharper 
and was shifted by 0.5 mm toward the center of the trap. 
In Fig. [H (b), the measurement at Kf = 200 V shown in 
Fig. [3] (b), obtained without filling the trap, is superim- 
posed for comparison. We note that, although the total 
number of ions in Fig. [3] (b) and that in the trap filled at 
Vi-f = 95 V are approximately the same, the distributions 
are different. 
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FIG. 2: Two-dimensional distribution of the Agj ion density. 
The total number A^'o of stored ions was about 5.5 x 10*. The 
eight circles represent the pole electrodes of the ion trap. 



In the above analyses of the absolute ion densities, 
the photodissociation cross section, a, was evaluated by 
Eq. ([2]) for each measurement. It was found to be a = 
(5 ± 1) X 10"^'^ cm^ at 415 nm for the present 300-K 
ion trap. The uncertainty represents a statistical error of 
the measurement; we mention that this evaluation may 
have an additional systematic error due to uncertainties 
in the estimation of the ion transmittance through the 
ion optics (see Section |IT1) and in the measurement of 
the laser pulse energy. Note that, as the spectral profile 
of photoabsorption is dependent on the temperature of 
the ions, the cross section at a given wavelength changes 
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FIG. 3: Radial distributions of the ion density obtained by 
one-dimensional scans on the horizontal axis for different 
amounts of loaded ions. The amplitude of the rf field Vi-f was 
200 V. Error bars show statistical standard errors of the mea- 
surements. Dashed lines show calculated distribution based 
on an adiabatic approximation according to Eq. (|10|) using 
the measured values for TVq. Upon convolution the calculated 
distribution with a Gaussian function of 2 mm FWHM rep- 
resenting the laser beam profile, we obtained the solid lines 
which should be compared with the experimental data. Ther- 
mal effects are not considered in these calculations. 



accordingly [36 



IV. MODEL CALCULATION 

To explain the radial distributions of the ion density 
dependent on the trapping conditions, we have performed 
a model calculation based on an adiabatic approximation 
[18| . Here, effects of the static potential produced by Lin 
and Lout are neglected because they are much smaller 
than the magnitude of the effective potential for the long 
trap we have used (40 cm). Thus we assume a cylindrical 
symmetry along the trap axis. The transverse motion of 
the ions in a multipole ion trap is described as a motion 
in an effective electric potential [30| expressed by 



*off(?') = 



p 



9K? 






(3) 



where g, to, 2p, i7/27r, and rg denote the ion charge, the 
ion mass, the number of poles, the frequency of the rf 
field, and the inscribed radius of the ion trap, respectively 
[l8|. In general, the equilibrium ion density n{r) in such 
a potential at a given temperature T, is described by 



i{r) = no exp 



kRT 



$cff(r) + $sc(r) 



(4) 



where uq and fcs are a normalization and Boltzmann's 
constants, respectively, and ^sc{r) is the electric poten- 
tial due to the space charge of the ion cloud [20, [30| ■ In 
the general case, Eq. (j4]) is non-linear because $sc('') is 
related to the local ion density by Poisson's equation. 
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FIG. 4: Experimental (symbols) and calculated (line) radial 
distribution profiles of the ion density for various rf amplitude 
Kt. (a) The trap is fully loaded with iVo = 1.1 x lO'* and 
1.2 X lO'^ ions at Kf = 300 and 200 V, respectively, (b) The 
trap is fully loaded with A'^o = 1-4 x 10* ions at Kf = 95 V 
and partially loaded with 1.6 x 10* ions at Kt — 200 V (the 
latter are the same data as in Fig. EJb)). 



However, if one assumes the space charge effect to be 
neghgible, Eq. (U) simphfies to 



n{r) — no exp 



knT 



*cff(?') 



(6) 



This approxiraation was apphed to explain an ion distri- 
bution in a 22-pole trap containing about 10'^ ions [3l[. In 
the present study, we have about 5 to 6 orders of magni- 
tude more ions so that we load the trap until saturation. 
Therefore we cannot neglect the space charge effect, and 
as a result we rearrange Eq. ^ as 

-$,,(r) = ^^ln^ + <i>eff(r). (7) 

q no 

At zero temperature the space charge due to the ion dis- 
tribution must exactly counterbalance the effect of the 
external potential $eff- That is to say, neglecting the en- 
ergy due to the thermal motion of the ions allows one to 
recast Eq. (H)) as [STJ 



^csir) + $sc(r) = 0. 
Equations ([5]) and ([8|) thus lead to 



n{r) = ^ V2$,ff(r), 



(8) 



(9) 



which determines the radial profile of the ion density. We 
will see later on that this approximation is not entirely 
fulfilled, but it is nevertheless sufficiently good to explain 
most of our experimental data. An alternative approach 
to the general derivation of the charge distribution in a 
cylindrically symmetrical external potential is described 
in the Appendix. Using Eq. ([3]) for the effective potential, 
we obtain the following radial distribution from Eq. ^: 



n{r) — p^{p 



^\2 ^oKf 
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2p-4 



(10) 



At T = we expect a sharp cut in the distribution at 
''max, which is determined by the number No of stored 
ions. In an octopole ion trap {p = 4), the ion density 
is proportional to r"'. Note that instead, in the case of 
a quadrupole (p = 2), Eq. ([T0|) predicts a constant ion 
distribution. In order to compare these predictions with 
the experimental data, the distribution n{r) was evalu- 
ated in the interval [0,rniax] and then convoluted with a 
Gaussian function, 2 mm FWHM, representing the laser- 
beam diameter. The results of the calculations are shown 
in Figs.[3]and|3]together with the experimental data. The 
maximal radius r^ax is such that the integral of n(r) over 
a 40-cm long cylinder of radius Tmax corresponds to the 
measured A^o- 

The standard way to create multipole potentials is by 
approximating the hyperbolic equipotential surfaces with 
cylindrical electrodes. In general, it is always possible to 
approximate a lower order multipole with the electrodes 
configuration designed for a higher order one. Thus, we 



approximated a quadrupole by wiring the neighboring 
poles of our octopole trap together to form four pairs 
of electrodes. The ion distribution measured for this 
quasi-quadrupole ion trap is shown in Fig. [S] The ions 
are concentrated around the axis with a rather flat den- 
sity profile as predicted by the present model calcula- 
tion. As the surface of the eight rods cannot perfectly 
follow the equipotential lines of a quadrupole potential, 
this quadrupole trap is not perfect and provides a smaller 
trapping volume than an ordinary four rods configura- 
tion. 
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FIG. 5: A two-dimensional distribution of ion densities mea- 
sured for AgJ ions in a quasi-quadrupole ion trap, which 
should be compared with Fig. [21 In this measurement the 
neighboring poles were wired pairwise together to form four 
pairs of electrodes. The total number of ions stored was about 
1.3 X 10* in this particular condition with Vrf ~ 360 V. The 
eight circles represent the pole electrodes of the ion trap, along 
with the polarity of the RF voltage applied to each electrodes. 



V. DISCUSSION 

The distributions calculated by using Eq. (ITU)) are com- 
pared with experimental results in Fig. [3] Note that three 
curves calculated for different Nq have the same profile 
except for rmax, which was determined to be 4.3, 3.0, 
and 2.4 mm for No = 1.2 x 10^ 1.6 x 10^, and 4.0 x 10^, 
respectively. The calculated curves reproduce well the 
overall features of the experimental results, particularly 
for larger A^o- This agreement indicates that the ring pro- 
files of the ion distribution are mainly governed by the 
space charge effect, which forces the stored ions toward 
the outer region. Note that the largest discrepancies be- 
tween the model and the experiment are found near the 
center of the trap and for a low ion density. This is due to 
the T = approximation, which is a poor approximation 
where f&eff and $sc are comparable to the energy of the 
thermal motion, i.e., at small r and/or small Nq. The 
room temperature, kgT = 26 meV, is equal to <I'cff('') at 
r = 2.2 mm for K-f = 200 V. 

Figures [21 and [S] show that the maximal ion densities 
in an octopole and in a quadrupole are similar, but the 



latter provides the maximum around the center of the 
trap, whereas in the former the maximal density is found 
on a ring at larger radius. Since it is easier to overlap a 
laser beam to an ion cloud that is concentrated around 
the center of the trap than to a ring, a quadrupole is 
probably the most favorable rf linear trap for laser spec- 
troscopy. However, since the phase space acceptance of 
a quadrupole is lower than that of every other higher or- 
der linear rf trap, it is possible that many ions injected 
into a quadrupole are lost. This makes the loading time 
of a quadrupole longer than for a higher oder multipole. 
In order to optimize the loading time and the laser beam 
overlap with the ions, one can imagine using a high-order 
multipole trap for loading and thermalizing the ions and 
then switching it to a quadrupole for spectroscopy. This 
can be done by controlling the potential on each electrode 
independently. 

According to Eq. ([TU|) . the maximum radius rmax 
within which ions are stably stored is determined by the 
total number of ions and by Vrf — for a given trap geom- 
etry and a given Jl. For full traps, however, two other 
factors reduce the value of rmax- One is the physical ge- 
ometry of the traps: according to Gerlich [l^ a typical 



geometric limit 



is 0.8 ro due to the space necessary 



for micro-motion wiggling. The other is the breakdown 
of the adiabatic approximation due to large rf amplitude. 
The adiabaticity is quantified by means of a dimension- 
less adiabaticity parameter r?(r), as described in Refs. 
fli and M: 



r]{r) =2p{p- !)■ 



qVri 



ro 



p-2 



(11) 



There is a maximum value of r/ = jymax for which the 
rf heating makes the motion of the ions unstable. Since 
r](r) increases with r, there is a critical radius Tc at which 
^(^c) = ??max- Whereas Tc decreases with increasing Kf, 
?'gcom is independent of Vif. Therefore, rmax is equal to 
TgcoTa at a sufficiently low V^f. As Vrf grows, the number 
of stored ions increases according to Eq. ([TU| until r^ is 



reduced to r„ 



A further increase of Vi-f causes a loss 



of the ions because rmax (now equal to Tc) is reduced. 
We summarize these two effects as r^ax = rninfrc, r„coml- 



This behavior has been observed by Mikosch et al. [23 1 
in their measurements of the trapping potential depth. 

The reduction of the maximum radius fmax observed 
when Vrf is changed from 200 to 300 V as shown in Fig.|3] 
is attributed to the decrease of Tc . For the V^f — 300 V 
case, ?7max is determined to be 0.13, where the value 
^max = 3.7 mm is obtained from Eq. (jlOp and from the 
measured Nq- We note that this value for rmax is in 
good agreement with the data shown in Fig. |4l Other 
groups found Tymax — 0.36 ± 0.02 for a 22-pole trap \23\ 
and 7?max ~ 0.2 for a 3D octopolc trap [30|- In addition, 
Gerlich suggested Tymax = 0.3 from numerical simulation 
|18| . Our value is much lower than the one for the 22-pole 
trap and the calculated one. The 22-pole experiment was 
done with less than 10"^ ions. Gerlich performed the sim- 
ulation for a single ion free from perturbation by other 



ions. In contrast, the present experiment was performed 
under a strong space-charge effect. The reduction of rymax 
indicates that the ion-ion interaction introduced an ad- 
ditional source of instability. 

On the premise that the same lymax = 0.13 be applied 
for the lower Vif, the values of re were calculated to be 4.4 
and 6.4 mm for V^i = 200 and 95 V, respectively. Clearly, 
rmax is limited by rgoom at T^f = 95 V. At T^f — 200 V 
we were able to store the maximum number of ions and 



therefore rr, 



the Tc value of 4.4 mm is 



in good agreement with the rmax value extracted from 
Fig. m The value of rgeom/ro = 4.4/5.5 = 0.8 is consis- 
tent with Gerlich's estimation of rgcom/ro < 0.8. 

The values of rgeom and Tymax are specific of this trap ge- 
ometry and possibly valid only in the high density limit, 
but otherwise independent of any other experimental pa- 
rameter. We have extracted these values using the mea- 
sured number of ions A'o together with Eq. (1101) . which 
describes the ion distribution in the high density limit. 
As the measurement of absolute numbers is always chal- 
lenging, it is interesting to consider the inverse problem: 
namely, the determination of the absolute number of ions 
in the trap, iVg, based on a relative measurement of the 
ion density and on the knowledge of the characteristics 
of the trap. This is only possible when the densities are 
high enough for the interaction between ions to become 
relevant in shaping the ion distribution. Then, one can 
turn a relative measurement of the ion distribution into 
a measurement of the interaction strength between ions 
and, thus, into a measurement of the absolute number 
of ions. Close to the space charge limit, one obtains Nq 
directly by integration of Eq. (|10p between and rj^^xj 
where r^J^^x is th^ maximum of the ion distribution, which 
is determined experimentally. 



VI. SUMMARY 

We studied the radial distributions of ions stored in 
a linear octopole ion trap near the space charge limit 
by monitoring photofragmentation yields as a function 
of the laser position. For the highest densities, we ob- 
served that the ion distribution has a ring profile. We 
showed that this is a typical feature of a multipole ion 
trap. The quadrupole potential, however, is an excep- 
tion in the family of the linear multipole rf traps as it 
induces a uniform ion distribution even when the space 
charge limit is reached. These observations are predicted 
and explained by a simple model based on equilibrium 
between the effective potential and that produced by the 
charge of the ions. The only approximation we used is 
that the energy related to the thermal motion of the ions 
inside the trap is negligible in comparison with the ef- 
fective potential generated by the multipole. This ap- 
proximation is fulfilled except for the regions where the 
effective potential is very flat, which, however, contain 
very few ions when the trap is full. 

The maximum adiabaticity parameter, rymax, was esti- 



mated to be 0.13, when the trap is fuU. This value is lower 
than those found in other studies under low ion-density 
conditions. We tentatively attribute the reduced value 
of 77max to an additional source of instability induced by 
repulsive forces among the stored ions. 

Understanding the space charge effects allows to extact 
the absolute number of ions in the trap and their absolute 
density based on the relative distribution. This provides 
a way of measuring absolute numbers without knowledge 
of absorption cross sections nor detector efficiencies. 
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Appendix: Charge distribution in a cylindrically 
symmetrical external potential 

The derivation of Eq. is presented in an alterna- 
tive way for a general cylindrically symmetrical external 
potential. The charge distribution is derived under the 
assumption that it is translational invariant along the z- 
axis, where it extends from negative to positive infinity. 
Thus, all extensive quantities as charge or energy repre- 
sent a "per unit length" value in the following. 

Assume that a certain amount of charge per unit 
length, Qoj is allowed to distribute freely in a cylindri- 
cally symmetrical external potential ^oxt. In cylindrical 
coordinates r, z, and 0, 



$ext(r-,z,6') = $cxtW- 



(A.l) 



The charge will distribute radially with density p(r) as to 
minimize the total electrostatic energy E per unit length 
of Qo in 'I'oxt- The density, p(r), is to be calculated as 
follows for any given external potential ^^y±(r). 

It is convenient to introduce the cumulative charge 
Qir), which is the amount of charge per unit length 
within a cylinder of radius r around the z-axis. Q{t) 
and p{r) are related by 



Pir) = ^g'M, 



(A.2) 



where Q'(r) is the derivative of Q(r) with respect to r 
and Q(r) is defined for < r < oo, its boundary values 
are Q(0) =0 and Q{oo) = Qo- The total energy E per 
unit length of the charge distribution is composed of the 
internal energy E^ of p(r) (i.e., the electrostatic energy 



resulting from the repulsion between the volume elements 
of p(r)) and the external energy, i?e, of the charge dis- 
tribution p(r) in the external potential. The latter is 
found by volume-integration of the product p(r)^ ci^\_(r) 
and can be expressed with the help of Eq. (jA.2p as a one 
dimensional integral along the r-coordinate: 

^e= /pW$ext(r)dA= / Q'(r)$extWdr, (A.3) 

J A Jo 

where A represents the volume in the r-9 plane. To ac- 
cess Ei{r) it is convenient to derive the strength of the 
radial internal field E from the first Maxwell equation 
applied over the surface of an infinite cylinder of radius 
r along the z-axis: 



E^9^e. 



(A.4) 



From classical electrostatics it follows that Ei is given by 



E = '-^ I E.EdA= r9Mdr. 



47r£o?' ' 



(A.5) 



Combining Eqs. (jA.3[) and (IA.5I) . the total energy can be 
written as a single integral over r 



"'L {«'('■)*»" +3} 



dr. 



(A.6) 



Equation (jA.Sp is already in the canonical form for the 
calculus of variations. The distribution Q{r) that min- 
imizes E under the given boundary conditions can be 
found by solving the Euler-Lagrange equation 



d dL dL 



dr dQ' dQ 



(A.7) 



where L = {Q'(r)$oxt(?') + %^], the integrand of 
Eq. (|A.7|) . is the Lagrangian associated with the vari- 
ational problem. Equation (jA.7p leads to 



Q{r) = 27reor$;,t(r), 



(A., 



which can be used to finally derive p{r) with the help of 
Eq. 1X2]): 



P{r) 



*cxt(0 



$;;,(r) = V2$ext(r) (A.9) 



which is equivalent to Eq. ^ and is applicable to calcu- 
late p{r) for any given cylindrically symmetric potential 

^oxtW- 
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